Transverse angular momentum of photons 
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We develop the quantum theory of transverse angular momentum of light beams. The theory 
applies to paraxial and quasi-paraxial photon beams in vacuum, and reproduces the known results 
for classical beams when applied to coherent states of the field. Both the Poynting vector, alias the 
linear momentum, and the angular momentum quantum operators of a light beam are calculated 
including contributions from first-order transverse derivatives. This permits a correct description of 
the energy flow in the beam and the natural emergence of both the spin and the angular momentum 
of the photons. We show that for coUimated beams of light, orbital angular momentum operators 
do not satisfy the standard commutation rules. Finally, we discuss the application of our theory to 
some concrete cases. 
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I. INTRODUCTION 

The quantum theory of light assigns a longitudinal 
component J • k/|k| = Ha of spin angular momentum 
(SAM) to a photon of energy fko and momentum fik, 
where a = ±1 for a circularly polarized photon, and 
(7 = for a linearly polarized one. At optical frequen- 
cies, however, the representation of a photon by a single 
plane wave mode of sharp angular frequency a; and wave 
vector k is quite unrealistic. Rather, a bona fide optical 
photon should be described as a wave packet formed by 
the superposition of many (possibly infinite) plane waves 
of different frequencies and wave vectors. As a result 
of this superposition a complex spatial structure of the 
photon field may be generated and tailored in order to 
carry an orbital angular momentum (0AM). Such pos- 
sibility was envisaged in 1992 by Allen, Woerdman and 
coworkers [l| who showed that a Laguerre-Gauss beam 
of light propagating in the z direction with a wave 
front of the form {x + iy£/\i\y^\ possesses a z coinponent 
of orbital angular momentum of M per photon [3|, with 
£ G {0 ± 1, ±2, . . .}. This result boosted the interest of 
the physics community for light beams with angular mo- 
mentum which had found numerous applications ranging 
from quantum cryptography Q to the realization of EPR 
entangled systems @ (see, e.g., ^ and Q for recent sur- 
veys). 

Previous authors have presented classical [1, Q and 
quantum [H, UH treatments of light beams with spin 
and orbital angular momentum. In these studies the at- 
tention was mainly devoted to the longitudinal (namely 
parallel to the beam propagation direction) component 
of the angular momentum. This was probably due to the 
fact that the spin angular momentum of photons can only 
be defined along the direction of propagation. However, 
it was very recently noticed that transverse, as opposed 
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to longitudinal, components of angular momentum may 
be responsible for interesting phenomena as, e^., the so- 
called geometric spin Hall effect of light fT2','T3"l. A clas- 
sical theory of transverse optical angular momentum was 
developed in p^ . 

In this paper we present a quantum theory of trans- 
verse angular momentum of photons. We apply the exact 
quantization scheme for light beams described in [13, HI] , 
to the development of a rigorous theory of quasi-paraxial 
photon fields, namely fields represented by the Lax et 
al. power series expansion [l^ truncated at first-order 
terms. The zero-order terms in the Lax expansion are 
exact solutions of the paraxial wave equation. However, 
it was shown in [3 that the presence of these terms 
solely is not enough to guarantee a correct description of 
both the energy flow and the spin angular momentum in 
the beam. Thus, we have included first-order transverse 
derivatives in our description of the photon fields. In this 
manner we were able to build a self-consistent theory of 
transverse angular momentum which displays some non- 
trivial characteristics as, e.g., anomalous commutation 
relations between angular momentum operators. 

This paper is structured as follows: In Sec. II we 
shortly review the exact quantization scheme [13, [13 and 
apply it to the present scenario. Then, in Sec. HI we 
derive closed expressions for both the linear and angu- 
lar momentum operators of the photon fields. By using 
these results, we show in Sec. IV that such operators do 
not fulfill canonical commutation relations in the paraxial 
regime of propagation, and discuss these findings. Sub- 
sequently, in Sec. V, we study some specific states of the 
fields that illustrate the occurrence of transverse compo- 
nents of the angular momentum. Finally, in Sec. VI we 
summarize our results. 



II. QUANTIZATION OF THE FIELDS 

In this section we illustrate the quantization procedure 
for quasi-paraxial beams of light. 
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A. Exact field quantization expansion of the positive-frequency part of the electric 

field operator in the Coulomb gauge can be written as: 

In Ref. 15] it was demonstrated that for a light beam 
propagating in the forward z direction, the plane-wave 



E+(r,t) = « / / dfcojdA:. 



dw 



1/2 



exp [~iuj [t — z/c) 



X ^ e^(k^, w)a^(kj_,a;) exp ikj_ • x - iz— (l - (fcj_,a;)) 



(1) 



r 



where r = (a;,j/, z) is the position vector and x = 
its transverse part in the plane z — const. Moreover, 
= {kx,ky) is the transverse part of the wave vector 
k — {kx, ky, kz), with fc = |k| = uj/c, and 

^^^,ik^,c.)^^i-'^y\ 

where k± = |kj^| = (fc^ + kyY^^. Note that the function 
Kz{kj_,uj) > is strictly defined only for uj > k±c, and 
that only the parts of the field propagating in the positive 
z direction are included in Eq. ([!]). 

The operator ap(kj^, w) annihilates a photon with wave 
vector k and polarization e^(k_L,a;), and satisfies the 
canonical commutation rules 

[a^(k^,c.),at,(kl,w')] = 5^^,5{k^ - k'^)6{io - oj'). 

(3) 

The three unit vectors |ei(k_L, w), e2(kj^, cj), k/fc} form 
a right-handed Cartesian frame that can be obtained 
from the reference basis {x, y, z} via a rotation around 
the axis 

n = i X k/\z X k|, (4) 

by the angle 9 between k and z: 

6 — arcsin(fc^/fc) = arcsin(fc^c/aj). (5) 

If we represent such a rotation by means of the Rodrigues' 
formula [Toj via the matrix 

R{9,n) ^cxp{eE) ^ I + E sine + E^{1 - cose), (6) 

where / is the 3x3 identity matrix and E denotes the 
antisymmetric 3x3 matrix of elements 
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E,j =~^e,jk{n)k, i,jG {1,2,3}, (7) 

k=l 



{sijk is the completely antisymmetric Levi-Civita sym- 
bol), then we obtain 

ei(k_L,w) =R{e,h)x, (8a) 
e2{k^,uj) ^R{e,h)y, (8b) 

and k/k = R {e, fi) z. It is worth noting that from Eq. 
([5]) and the definition ([2]), it follows that 

sin^? = fc^c/a; ^^(fc^,^) = cos0. (9) 

The geometry of the "global" i} and the "lo- 

cal" {ei(kj^, w), e2(kx, oj), k/fc} Cartesian frames, is il- 
lustrated in Fig. 1. 



y 




FIG. 1: Illustrating the geometry of the problem. 



B. Quasi-parcLxial quantization 

Equation ^ is exact, therefore it can be used to de- 
scribe any field that propagates in the positive z direc- 
tion. However, most optical experiments use narrow- 
band and well-coUimated beams which satisfy, respec- 
tively, the conditions 

Aw<wo, and 6*0 < 1, (10) 
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where ujq is the central frequency of the bandwidth Aw, 
and 00 is the angular spread of the beam around the cen- 
tral wave vector ko = zluq / c. For states of the radiation 
field whose excitation bandwidths and angular apertures 
satisfy Eq. (|10p . the k-space in Eq. ^ can be restricted, 
without significant error, to the intersection I between 
the cone of axis ko oc £ and aperture 29q, and the spher- 
ical shell of radius wq/c and thickness Auj/c, as shown in 
Fig. 2. For any k G I we have do 1, and from 



of the electric field operator can be written as 





1 











'1 



FIG. 2: The grey area represents the domain X within the 
k-space. 



Eqs. ([5|9| it follows that 

K,{k^,Lu) ~ l-isin^e = l-i(fc^c/c^)^ (11) 

where k±c/uj <^ 1. Thus, the function Kzik±,uj) inside 
the square root in Eq. ([T]) can be approximated by 1, 
while the second order term must be retained in the z- 
dependent part of the last exponential in Eq. ([1]) to ob- 
tain a nonzero result: 



-iz-(l - iq,uj)) ~ ^. 

C ZU! 



(12) 



Moreover, we can approximate Eq. ([6|) with 

R{e,n) ~ I + EsinO = I + {k^c/u})E, (13) 
and from Eq. ([8|) it readily follows that 



ei(k^,w) ~ A - zk^/k, 
e2(kx, w) ~ y - zky/k. 



(14a) 
(14b) 



E+(r,t)=i/ duj i ^Q^s ) exp[-iuj{t- z/c)] 



r ^ 



X exp I ik^ • X — iz 



2^y ' 



(15) 



where d^fcj^ = dk^dky, xi = x, X2 = y, and we inter- 
changed the order of integration. Equation (fT5|) can be 
further simplified by noting that 

d 

-k^ exp {ik± ■ x) = I- — exp {ik± ■ x) , (16) 
dx^ 

which permits us to rewrite Eq. (llSp as 



E+(r,t) / duj 



V 47r£i 

2 



— I exp\~iuj{t — z/c)] 
qcJ 



X ^ (^A^ + • V_L^ ap(x, 2:,w), (17) 



where we have defined Vj_ — xd/dx + yd/dy, and 

a^{x,z,uj) ={2t:)~^ j d?k±_ a,,(k_L,w) 

X exp ikj^ • X — iz^^ ) . (18) 
V 2cjy 

At any plane z = const., this is still a bona fide quantum 
harmonic oscillator annihilation operator, as it satisfies 
the following commutation rules: 

[a^(x, z, i:j),aj^,(x',z,w')] = 5^^,5{y. - x!)6{uj - w'). 

The integral ^TEh can now be evaluated by using the fol- 
lowing relation [lli] 

-exp(^.k^.x-zz— 

^^'^nni{^±:^)i'n7n{x,Z,Uj), (19) 



where V'nm(x, z, w) is a complete orthonormal set of func- 
tions on x. 



V'*™(X, Z, LO) Ipn'm' (x, Z, Uj) d^X = 5 nn' Srayn' , 
^^lmi^,Z,Uj)^Jrim{y^ ,Z,Uj) =(5(x-x'), (20) 



Finally, for excitations of the field satisfying Eq. ([T0| . 
we can extend the integration over w in ([T]) from to cxd 
without relevant error, and the positive-frequency part 



and n,ni G Z are the appropriate integer labels for the 
set. In order to fulfill Eq. (|19p . the functions ?A„m(x, z, uj) 
must be chosen amongst either the Hermite-Gauss (HG) 
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or the Laguerre-Gauss (LG) sets of solutions of the parax- 
ial wave equation [2] , and 



V'nmlkijw) = (27r) ^ J ^Anmlx, 0,a;)exp(-ix- k^)d^a;, 

(21) 

is just the Fourier transform of '0„m(x, z, evalu- 
ated at 2; = 0. Using ([19]) inside (fT8|) we obtain 

a^(x, z,a;) = ^ a^„,„(a;)V'„m(x, z, w), (22) 

n,m 

where the operator 

a^mm(w) = / d^fcj_?A*„j(k^)a^,(ki,tj). (23) 



annihilates a photon with polarization in the spatial 
mode ipnm- From Eq. (jS]) and exploiting the orthogonal- 
ity of the modes ipnm , it is easy to verify that 

[a^mm(w),aJj/„/„/(w')] ^ 6^^'6nn'Smni'S{uJ - Uj'). (24) 

Finally, we can write the positive-frequency part of the 
electric field operator as 

E+(r,t) = i J doj yj:^^J exp[~iuj{t ~ z/c)] 

(25) 



It should be noticed in the expression above the presence 
of the transverse gradient V_l which is equivalent to the 
first-order term in the Lax et al. expansion [Tsj . 

By following the same route that lead us to (|25|) . it is 
not difficult to see that within the same approximations 
([T0| the positive-frequency part of the magnetic field op- 
erator can be written as 

i f°° f hijj \ ^^'^ 
B+{v,t) ^-j^ dio ( ) exp[-zw(t- z/c)] 

(26) 

where the symbol " x " denotes the ordinary vector cross 
product. 

C. Monochromatic limit 

In the laboratory practice, one often deals with laser 
beams whose bandwidth is so narrow that they can be 
considered basically monochromatic. For this case the 
formalism that we have developed above may be redun- 
dant and simplified expressions can be used. In order to 
pass from the general case above to the monochromatic 



limit it is convenient to make, as a preliminary step, the 
passage from a continuous to a discrete frequency spec- 
trum by letting 



UJj = LOo +jAuj, 



(27) 



where je{0,±l,±2,...}. In this limit the continuous- 
frequency annihilation operators are transformed to the 
discrete-frequency ones via the rule 



-1/2, 



(28) 



and integrals over continuous frequency are converted to 
sums over the discrete index j according to 



(29) 



The discrete-frequency electric field operator is obtained 
by applying the rules (1271129^ in ((25)) . thus obtaining 



E+(r,i) = *5] 



2eoc/Aa; 



1/2 



exp [—iu!j{t — z/c)] 



E 



X > a 

/i,n,m 



(30) 



where the discrete- frequency annihilation operator a^nmj 
satisfies the commutation rules 

[S/irmy, ajj'n'm'j'] = S f^fj.' ^nn' Smm' Sjj' ■ (31) 

For excitations of the radiation field that satisfy ^TU\\ we 
can keep the solely term j = in (j30| to obtain the strict 
monochromatic limit 



E+(r, t) = i 



2£oc/Aw 



1/2 



exp [—iu!o{t — z/c)] 



Opnm (Xf^ +iz—X^ ■ V± ] V'«m(x, Z, Wq), 

(32) 



^/inmj I Q • It 



where we used the shorthand ap„,„ 
is worth noting that in the expression above, the strict 
paraxial limit would be obtained only by neglecting the 
term proportional to a;^ • V^. 

The same procedure that lead to Eq. ([5^ can be fol- 
lowed to obtain the following expression for the positive- 
frequency part of the magnetic field operator in the 
monochromatic limit: 



B+(r,i) 



hujo 



2eoc/Au} 



1/2 



exp [~iu!Q{t — z/c)] 



p,nm I -2 X Xp^ + i X X V_L ) V'; 

Wo 



mn»,X, Z, Woj. 

(33) 



5 



III. LINEAR AND ANGULAR MOMENTUM OF 
THE FIELD 

In the previous section we have derived expUcit expres- 
sions for the electric and magnetic field operators. This 
allows us now to calculate both the linear and the angular 
momentum of the quantized electromagnetic field. 



A. Linear momentum 

The normal-order linear momentum density operator 
: P{r,t) : is equal to 1/c^ the Poynting vector normal- 
order operator : S(r,t) : and it is expressed in terms of 
the electric and magnetic field operators as [20j : 



:P(r,i): = 



E(r, t) X B(r, t) - B(r, t) x E(r, t) 



(34) 



(35) 
(36) 



If in the expression above we substitute 

E(r,t) =E+(r,t)+E-(r,i), 
B(r,i) =B+(r,i)+B-(r,0, 
with 

E-(r,i) = [E+(r,^)]^ B-{r,t) = [B+{r,t)]\ (37) 
we can write Eq. ([M)) as the sum of 8 terms: 



:P(r,i): 



{( 



E" X B+ - X E+ 
E+ X B" - B+ X E 



E" X B" 



E^ 



B X E 



B+ X E" 



-)}:, (38) 



where, for sake of clarity, we have omitted the explicit 
space and time dependence of the electric and magnetic 
field operators. The last four addenda of this sum con- 
tain integrands with terms as a^nm(w)a^'„'m'(w') and 

'^lmmi'^)^''fj.'n'm'(^') ^hich osciUatcs at frequencies higher 
than iM : ujQ — Auj < uj < ujq + Aw where, once again, ujq 
is the central frequency of the beam and Aw its band- 
width. As most optical detectors integrate the received 
signal over a time interval T much longer than l/Aco, the 
last four terms in Eq. p8p can be neglected in the limit 
T ^ l/Aw. Within this assumption, the time-integrated 
power flow trough the detector surface (coincident with 
the x, y plane) or, equivalently, the linear momentum op- 
erator per unit length, can be written as: 




P(r,i)dt \dxdy: , 



(39) 



where the spatial integration is extended over the whole 
plane. Although the integration time T in Eq. ([M)) 



is finite, the integration interval [— T/2,T/2] can be for- 
mally estended to (— c», cxd) without significant error for 
narrow-band beams. Thus, substitution in Eq. (pQ)) from 
Eqs. (EMMl) gives 

p^_/ dwfiwK:^:^:^""-'™' 

^ ^u— 1 n,m n'm' 



X aLm(w)a^n'm'('^) k 



(40) 



where 



_2. 

2i 



ip*^ (z - iXV^) -ipn'm' dxdy 



-'n' ,n+l 



nn' ^mm' ■> 



(41) 

with A — 2ttc/uj. Not surprisingly, Eq. (j40|) has the same 
form of Eq. (21) in Ref. [l^. The main difference is in 
the form of the quasi-paraxial linear momentum oper- 
ator z — iX^ ± as compared with the truly "quantum- 
mechanical" linear momentum operator —iK^ in coor- 
dinate representation. This difference is substantial and 
will manifest its effects later, when we will calculate or- 
bital angular momentum operators and their commuta- 
tion relations. The strict paraxial limit is obtained by 
tacking the limit ^ in Eq. PO]) . In this case only 
the longitudinal component Vz keeps a nonzero value, 
but this is inconsistent with a correct representation of 
the energy flow in the beam [18]. 

Finally, the three components of the operator P can 
be calculated explicitly by substituting Eq. (|4T|) into Eq. 
(|40l) obtaining 



(42) 



X 



aLm('^)aM":™+l(^) - aL,m+i('^)aM"m(w) h (43) 



for the transverse components, and 

•^•^ fji,= ln,m 

for the longitudinal one, which is simply proportional to 
the total number of photons per unit length of the beam. 
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B. Angular momentum 



Again, substitution in Eq. (|46|) from Eqs. (|25II26I) gives 



The calculation of the normal-order angular momen- 
tum density operator : 3{r,t) : proceeds along the line 
delineated in the previous subsection, starting form the 
standard definition 



:J(r,t): = :rxP(r,0:, 



(45) 



and arriving to the time-integrated angular momentum 
operator per unit length: 



X aLm(^^)aM'"'m'(w) 



finm, fi n m 



(47) 



1 



T/2 



T/2 



3{r,t)dt dxdy 



(46) 



where we have defined 



j^nm,^'n'm' ^ j ^.^^^^ X + V r X (z - iXV ^n'm' dxdy 



{nm'Y/'^5n,n' + l57n' ,m+l " ("-'w)^/^5„/,„+i(5m,„/ + i 



(48) 



with X ic^/ = = zct^^/. From Eqs. (|41l48p it 

follows that we can write 



J = S®Il 



(49) 



where X5 and Xl are the identity operators in the spin 
and orbital angular momentum spaces, respectively, and 
we have defined 



(5.Z) 



£ = r X -P, 



(50) 
(51) 



with e {1,2}, and «S • x = = <S • y. This result 
is the quasi-paraxial analogous of Eq. (20) of Ref. [l^, 
and it shows the separation of the total angular momen- 
tum of the beam in its spin and orbital parts. In the 
paraxial limit ^ 0, Eq. (|^5| is dominated by the 
transverse part and both the spin and the angular con- 
tribution to Jz become negligible. This is in agreement 
with the results of Haus and Pan [l^ who have shown 
that a self-consistent description of angular momentum 
of light beams cannot be achieved in a purely paraxial 
context. 

From Eq. (pS)) the three components of the angular 
momentum operator per unit length can be explicitly cal- 



culated obtaining 
1 A 



Jx — 



0-^0 U = l„.m 



/2 



aLrn('^)aM":™+i('^) + aL,m+i(^)aM"™('^) T' (52) 



J. 



^ ^ u—l n.m 



(53) 



and Jz = Sz + Lz, where we have defined the spin and 
the orbital angular momentum longitudinal components, 
respectively as 

,2^7 /_ dtj/ia;^ [^a|„„(w)a2„m(w) 



Sz 



c^T i Jo 

-4nm('^)«l«™('^)] ' 



(54) 



and 



1 A 



c2T 



A f 

- dwfiw V [{n + l){m + l)] 



1/2 



fi.n,m 

(55) 
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IV. COMMUTATION RULES 

At this point we have collected all the ingredients nec- 
essary to calculate the commutation relations between 
the linear and angular momentum operator components. 
From Eqs. (|40l47p it follows that both these operators 
have the form: 



Ua 



1 

— dcohwJ2&li^a^'aA', (56) 

-^0 A,A' 



where U e {J, P}, U e {J,V}, and a e {x,y,z}. In ad- 
dition we have introduced the cumulative labels A; A'] . . . 
that embody the three indices /i, n, m; /i', n', m'; . . . re- 
spectively, so that, e.g., J2a = E^E„E„- Now, we 
assume the validity of the following commutation rela- 
tions: 



(57) 



where U,V,W S {J,P}, U,V,W S {J,V}, a,b,c G 
{x,y,z}, and are numerical coefficient to be deter- 
mined and summation over repeated indices is under- 
stood. Then, it is not difficult to see that 

^ A,A' 



2ll fa 

cT 



A,A' J 



^Jabc c 
CI 



(58) 



where the factor cT in the last line of the equation above, 
plays the role of "natural" unit length for the problem 
under consideration. We remind that the validity of Eq. 
([55]) is subject to the assumption ([57]) that must be still 
verified. We begin such check by noticing that from Eq. 
(HTI) it is easy to calculate 



\Pa-,Vi\ =0, a, 6 G z}, 



(59) 



which has the form (j57p with = 0. Similarly, from 
Eqs. (|48I49P it follows that 

[Ja, J&] ^T® {U.tlh\ + [5a,5fo] ®X. (60) 

which reduces to its first term solely sinc6 Sx — — Sy 
trivially implies 

[5a, 5^] =0, a,h^{x,y,z). (61) 

By using Eq. ([i5|), an explicit calculation furnishes 



(62) 
(63) 
(64) 



that, together with Eq. ^ gives /f^f = ea6c(l - Scz), 
that amounts to a violation of canonical commutation 
relations for the total angular momentum. Before dis- 
cussing this somewhat surprising results, let us conclude 
the calculations by showing that 



\_^x-> 'Py 
\_^y 1 Px 



-0, 
= 0, 

= iXVy 



(65) 
(66) 
(67) 
(68) 
(69) 
(70) 



which follows from Eqs. (|4I|48M9p and from the triv- 
ial identity [7^a,56] = 0, as P and S operates upon 
different linear spaces. For short, Eqs. (I65ll70p give 
/fhc = £a6c(l - Sb 



A. Discussion 

In order to discuss the results present above, it is use- 
ful to adopt the notation O to indicate the operator O 
in non-relativistic two-dimensional quantum mechanics 
in the Schrodinger picture. Thus, for example, x = x, 
p — —ih'V±, etc. Note that in this representation the 
longitudinal coordinate z is not a dynamical variable, 
but a parameter that plays the role of "time" in the 
Schrodinger equation [21|. If with \ti), ... we denote 
the eigenstates of a one-dimensional harmonic oscillator, 
then from Eq. (|4T|) it follows that 



Vx oc {n\px\n'){m\m'), 
Vy oc {n\n'){m\py\m'), 
Vz oc {n\n'){m\m'). 



(71) 
(72) 
(73) 



The equations above show that while Vx and Vy behaves 
like Cartesian components of the canonical linear mo- 
mentum, Vz does not as it is proportional to the identity 
operator. In a similar manner we can see that Eq. (|48p 
implies 

Cx cc {n\y\n'){m\m'), (74) 
Cy OC — {n\n'){m\x\m'), (75) 
Cz (x{n,m\xpy ~ypx\n',m'), (76) 

namely Cx ^ y, Cy ^ ^x and Cz ^ L^, where the sym- 
bol stands for "behaves like" , and Lz — xpy — ypx ■ 
If we look at Eqs. (|7m76l) the origin of the anomalous 
commutation relations (I62II70|) becomes clear, since 

[Cx,Cy] ^ [y,~x] =0, (77) 
[Cx,Cz] ^ [y,Lz] ^ ihx ^ -iXCy, (78) 
[Cx^Vy] ~ [y,py\ =ih^ iXVz, (79) 
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and so on. In physical terms, the essence is that any 
well collimated beam is basically an eigenstate of the z- 
component of the linear momentum operator which, for 
such beams, practically reduces to a c-number. Thus, 
in order to recover canonical commutation relations, it 
is necessary to deal with beams with either high angu- 
lar aperture for whose our first-order approximation 
breaks down, or with a direction of propagation that de- 
viates from the reference axis z by an angle grater than 
^0 [13 ■ a final remark, it should be noticed that viola- 
tions of canonical commutation relations for the angular 
momentum of an electromagnetic field of arbitrary shape, 
were already reported by van Enk and Nienhuis [lo| . 



V. EXAMPLES 

In the previous section we have completed the study of 
the formal properties of the linear and angular momen- 
tum operators for quasi-par axial beams. In this section 
we will illustrate the usefulness of our treatment by ap- 
plying it to the common case of coherent excitations of 
the electromagnetic field. 

Let us begin by considering the quantum coherent state 
representing an ordinary monochromatic laser beam pre- 
pared on the laboratory bench in the spatial mode 
V'nm(x, Wo) and linearly polarized along the direction x^: 



exp I aa 



it 



a aA 



|0), 



(80) 



where A = fj..,n, m. For non-monochromatic coherent 
states of the field, this simple expression immediately 
generalize to [l^l 



[aA]) =D[aA]\0), 



(81) 



where 

D[aA] = exp 



do; 



(^aA{uj)a\{uj) - a*A{uj)aA{i^ 



(82) 



is the displacement operator for the multi-mode quasi- 
paraxial electromagnetic field, and square brackets [ua] 
indicate functional dependence. It is easy to see via a 
direct calculation that the amplitude functions aA{i-o) are 
the eigenvalues of the annihilation operator 0^1(0;): 



aA{uj)\[aA']) ^ aA{iL>)SAA'\[aA]), 



(83) 



where Saa' = ^ fj.fj.' Snn'S mm' ■ The coherent state \[o:a]) is 
an eigenstate of the positive- frequency part of the electric 
field operator [23j 



E+(r,i)|[aA]) =E+[a^](r,i)|[aA]), 



(84) 



where the eigenvalue E^[aA](r, t) is the analytic signal 
[20} of the classical field generated by the excitation with 



spectral amplitude a a = a^nmi^) of the mode ipnm- 

/•oo / \ 1/2 

E^[Q!A](r, t) = i / do; I ) exp[-iuj{t - z/c)] 



X a 



(85) 



Multi-mode field coherent excitations are handled exactly 
in the same manner by writing the iV-mode coherent 
state l^fAr) = as 



N 



= \[aAMo^A.\) ■■■■■ = n (86) 

i=l 

and Eq. (|84l) becomes 



E+(r,t)|vI/^) =E+[W}](r,t)|vI/ 



N1 



(87) 



where 

N 

= ^E+KJ(r,t) 



i=i 

/•oo / \ 1/2 

i I Au [ I exp[— ici;(t — z/c)] 

7o V47reocy ^ ' ' 



fi^n,m 



Equations (|86ll88p describe the more general coher- 
ent state excitation for a quasi-paraxial beam. The ex- 
pectation value of linear and angular momentum opera- 
tors with respect to these states can be written, after a 
straightforward calculation, as 



l.N 



1 

dujhuj{d{uj)Ma{^))j,, (89) 

where U G {J,P}, U £ {J^V}, and we have used 
the suggestive notation (a(a;), /3(w))^ to indicate the 
scalar product between the two TV-dimensional vectors 
d{uj) = (a^i, . . . ,aA„) and f3{uj) = (/J^j, . . . , /3^„). For 
a monochromatic beam of central frequency uiq-, Eq. (j28p 
requires aA(w) — {Au;)~'^/^a^{u}o) (the superscript D 
stands for "discrete"), and Eq. ()89|) reduces to 



(90) 



It is instructive to apply Eq. ([90]) to a concrete case 
in order to see the physical meaning of the equations 
above. With this aim, let us consider a spatial mode of 
the field of the form u/(x, z, ujq), where u = ^x-f Tyy, with 
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|^P + |77p = l, isa unit vector that fixes the polarization 
of the beam, and 

/(x, z, ujq) = aooV'oo + aioV'io + aoiV'oi, (91) 

where xj^nm — ipnmi^, z, Wo) indicates the (n, m) Hermite- 
Gaussian mode [2l, and jaooP + l^ioP + |aoiP — 1 be- 
cause of normalization. The corresponding coherent ex- 
citation I^'at^g) can be written as 

1*6) =|[aioo])|[aiio])|[aioi]>|[a20o])|[a2io])|[a20i]) (92) 
= iCaoo) IC"io) l^aoi) haoo) \mw) 



where we have defined ai„m(wo) = a,, 
Eq. dMl) from Eq. (gl]) gives 



Substitution in 



(^-elPI^-e) = ^(iP. + yPy + iP.), (93) 



with 



'o Im (aSo"io) , 

00*^ 



Py =6'oIm(aooQ;oi) , 



P. 



1. 



(94) 
(95) 
(96) 



The expressions above are rich of information. First, Eq. 
([M)) shows that P^ = |q!ooP -I- \aw\^ + |aoiP = 1 fur- 
nishes the total intensity of the beam, and that it is 
unbiased with respect to the modes ijjnm, thus reveal- 
ing both its "identity" character. Second, if we define 
tan 6*^ = P^/P^, (// e {x,y}), then from Eqs. (|Mlig5)) 
it follows that the spatial mode /(x, z, wq) deviates from 
the axis z only when ago 7^ 0, namely only when the fun- 
damental Gaussian mode V'oo is present in such super- 
position. In other words, in order to define an "absolute 
deviation" it is necessary the presence of a reference mode 
V'nn which is symmetric with respect to the simultane- 
ous inversion x— y— j/. Finally, when aoo S 
Eqs. (|94ll95p reproduce the well-known result that either 
aio or aoi must have an imaginary part to guarantee a 
nonzero tilting angle 9^ {23]. 

In a similar manner we can evaluate the expectation 
value of the angular momentum operator obtaining 



(*6|J|*6 



^0 



[xJ^ + yjy + zJ^ 



where 



Jx = Wo Re (aooftoi) , 
Jy ^ - wq Re (aoo"io) , 
J^^X[a -I- 2Im(aioaoi)] 



(97) 



(98) 
(99) 
(100) 



with a = i{(,rj* — (,*r]) denoting the elicity of the beam. 
Here, if aoo S M, Eqs. (|98l) and (|99| imply that in order 
to have a nonzero transverse orbital angular momentum 
either aio or aoi must have a real part. However, it 
is well known that a superposition with real coefficients 
of the fundamental mode "000 with either ipiQ or de- 
scribes approximatively a displaced Gaussian beam along 



the X- or the y-axis, respectively [2J]. In other words, 
a lateral displacement of a Gaussian beam changes its 
transverse angular momentum or, vice versa, the occur- 
rence of non zero transverse components of the angu- 
lar momentum cause a transverse displacement of the 
beam [l2| . On the other hand, it is also known that 
a transverse displacement cannot affect the longitudinal 
angular momentum [25. 26] , as it is confirmed by Eq. 
poop which goes to zero when both aio and aoi are real 
numbers. However, for a pure Laguerre-Gaussian beam 
we have /(x,z,a;o) = V'fcii^p^o = (V'lo ± #oi)/\/2, 
namely aoo — ^ Jx = = Jy Thus, in this 
case 2Im(a*oaoi) = ±1 and Eq. (|100p furnishes Jz — 
X{a ± 1) which, in agreement with previous calcula- 
tions [1], shows that a Laguerre-Gaussian beam possesses 
1^1 = 1 units of orbital angular momentum along of the 
direction of propagation. 

It worth noting that the results of this section for a 
coherent beam are in perfect agreement with the classical 
results presented in Ref. pL^ . 



VI. SUMMARY 

In this paper we have applied the theory of quantized 
light beams to investigate the properties of the trans- 
verse components of the angular momentum operator of 
the electromagnetic field. It is known that for either mas- 
sive particles or photons localized in wave packets it is 
meaningful to talk about the total angular momentum of 
the system under consideration. This is evaluated as the 
integral over the whole K'^ space, of the angular momen- 
tum density of either the particle or of the field. However, 
when dealing with photons in beam-like states, the rele- 
vant quantity which can be actually measured on a labo- 
ratory bench, is the angular momentum per unit length. 
This is evaluated, in a plane of equation z = const, per- 
pendicular to the main direction of propagation z of the 
beam, as the integral over the transverse xy-plane of the 
angular momentum density of the field. Since in this case 
the integration does not extend over all 3D space, the an- 
gular momentum per unit length is not independent of 
time Q. Therefore, it becomes necessary to average this 
quantity over the measurement time T, and one is led 
to the expression shown in Eq. We have explicitly 

evaluated this time-averaged angular momentum per unit 
length J within the framework of paraxial optics, but in- 
cluding contributions from first-order transverse deriva- 
tives of the electric and magnetic fields. This inclusion 
permitted us to achieve a self-consistent description of 
both the energy fiow (linear momentum or Foynting vec- 
tor) and the spin and orbital angular momentum of the 
beam which appear to be naturally separated in the sum 
J = S + L. Then, we have calculated the commutation 
relations between the Cartesian components of J, L, S 
and we have found that they differ from the standard 
one. In particular, while is still a bona fide generator 
of rotations around the propagation axis z, the transverse 
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components Lx and Ly commute and, as it was already 
found at a classical level [T2| , they are strictly connected 
with the transverse coordinates y and —x of center of the 
beam, respectively. Finally, as a realistic example illus- 
trating the above mentioned connection, we calculated 
the expectation value of J between multi-mode coherent 
states of the electromagnetic field. It is well known that 
for these states quantum and classical expectation values 
basically coincide. Indeed, we found full consistency be- 
tween our results and the classical ones [13], namely we 
found that {L^) oc Aj, and {Ly) cx — A^,, where A^, and 
Aj^ are the transverse displacements of the center of the 



beam with respect to the propagation axis z. Further in- 
vestigations of the relations between angular momentum 
and beam shifts in classical optics, are illustrated in Refs. 
I27I l28i and references therein. 
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